A Super-Earth Orbiting the Nearby Sun-like Star HD 1461 

Eugenio J. Rivera 1 , R. Paul Butler 2 , Steven S. Vogt 1 , Gregory Laughlin 1 , Gregory W. 

Henry 3 , Stefano Meschiari 1 

Received ; accepted 

draft: hdl461_rev3.6 



1 UCO/Lick Observatory, University of California, Santa Cruz, CA 95064, USA 
2 Department of Terrestrial Magnetism, Carnegie Institution of Washington, 5241 Broad 

Branch Road, NW, Washington, DC 20015-1305, USA 

3 Center of Excellence in Information Systems, Tennessee State University, Nashville, TN 

37209, USA 



- 2 - 



ABSTRACT 



We present precision radial velocity (RV) data that reveal a Super-Earth mass 
planet and two probable additional planets orbiting the bright nearby GOV star 
HD 1461. Our 12.8 years of Keck High Resolution Echelle Spectrometer pre- 
cision RVs indicate the presence of a 7.4M ffi planet on a 5.77-day orbit. The 
data also suggest, but cannot yet confirm, the presence of outer planets on low- 
eccentricity orbits with periods of 446.1 and 5017 days, and projected masses 
(Msini ) of 27.9 and 87.1 M®, respectively. Test integrations of systems consis- 
tent with the RV data suggest that the configuration is dynamically stable. We 
present a 12. 2- year time series of photometric observations of HD 1461, which 
comprise 799 individual measurements, and indicate that it has excellent long- 
term photometric stability. However, there are small amplitude variations with 
periods comparable to those of the suspected second and third signals in the RVs 
near 5000 and 446 days, thus casting some suspicion on those periodicities as 
Keplerian signals. If the 5.77-day companion has a Neptune-like composition, 
then its expected transit depth is of order d ~0.5 mmag. The geometric a priori 
probability of transits is ~8%. Phase folding of the ground-based photometry 
shows no indication that transits of the 5.77-day companion are occurring, but 
high-precision follow-up of HD 1461 during upcoming transit phase windows will 
be required to definitively rule out or confirm transits. This new system joins 
a growing list of solar-type stars in the immediate galactic neighborhood that 
are accompanied by at least one Neptune (or lower) mass planets having orbital 
periods of 50 days or less. 

Subject headings: planetary systems - stars: individual (HD 1461) 
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1. Introduction 

Over 400 extrasolar planets are now known. The majority of these have been discovered 
by using precision radial velocities (RVs) to detect the reflex barycentric motion of the host 
star. We have had a large sample of over 1000 nearby stars under precision RV survey for 
the past 13 years at Keck with the High Resolution Echelle Spectrometer (HIRES). One o 



the t arget stars is HD 1461, a nearby GOV star, only 23.4 ±0.5 pc away (jPerryman et al. 



19971 ). This star has been on the Keck program since 1996 October. Over the past 12.8 
years, we have accumulated a total of 164 precision RVs that indicate a system of at least 
two planets orbiting this star. In this paper, we present all of these RV data and discuss 
the planetary system that they imply. 



2. Basic properties of the host star HD 1461 



HD 1461 (HIP 1499, HR 72, GJ 16.1, GJ 9009, BD -0838, SAO 128690, SPOCS 13) is 
a bright (V=6.46) and well-studied star. It h as been characterized in a n umber of stu dies, 



i nclud ing those of 



Valenti fc Fischer 



(120051 ) 



Takeda et al. 



(120071 ) . and 



Sousa et al. 



( 120081 ) . Table [T] summarizes recent determinations of the fundamental stellar parameters for 



HD 1461. Taken together, these properties indicate that it is an old, metal-rich, inactive star 
well suited for precision RV plan et searches. Using Ca H+K measurements taken between 



1994 and 2006. 



Hall et al. 



(120071 ) found that HD 1461 is one of the 13 targets for which the 
observed variability is zero within th e uncertaint i es. In 51 observations over 8 seasons, they 
measured a mean log-R^K °f — 5.04. 



Hall et al. 



(120091 ) found a mea n value of —5. ove r 



Wright et al. 



(12004 ) 



seven seasons with seasonal averages ranging from —4.96 to —5.01. 
found a mean value of —5.03, along with an estimated slow rotational period of 2 9 days 
and a n age of 6.3 Gyr. Our measurement of log -R^ K =— 5.00 leads to an estimate (jWright 



20051 ) of 1.59 ms 1 for the expected RV jitter due to stellar surface activity. The age of 
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HD 1461 was e stima ted as 4.2+2 x Gyr by 



Valenti fc Fischer 



(120051 ) and 7.12+};^° Gyr by 



Takeda et al. 



(120071 ). The large discrepancy and uncertainties for these chromospheric ages 



is not surprising since the correlation bet ween age and chromosphe ric activity becomes very 



weak for ages bigger than 2 Gyr (see e.g. 



Pace fc Pasquini 



((200J)). In summary, HD 1461 



is a nearby, bright star with physical properties that are quite similar to our own Sun, and 
is an ideal candidate star for the application of high-precision Doppler- velocity monitoring. 



3. Observations 



The HIRES (jVogt et al. Ill994l ) of the Keck I telescope was used to monitor HD 1461. 
A total of 164 Keck observations were obtained, from 1996 October 10 to 2009 August 10, 
a data span of 4687 days. The median internal velocity uncertainty for these Keck data is 
0.72 ms" 1 . 



19961 ) by placing an 



Doppler shifts were measured in the usual manner (IButler et al. 
iodine absorption cell just ahead of the spectrometer slit in the converging beam from the 
telescope. This gaseous iodine absorption cell superimposes a rich forest of iodine lines on 
the stellar spectrum, providing a wavelength calibration and proxy for the point-spread 
function (PSF) of the spectrometer. The iodine cell is sealed and temperature-controlled to 
50 ± 0.1 °C such that the column density of iodine remains constant. For the Keck planet 
search program, we operate the HIRES at a spectral resolving power R m 70, 000 and a 
wavelength range of 3700-8000 A, though only the region 5000-6200 A (with iodine lines) is 
used in the present Doppler analysis. The iodine region is divided into ~700 chunks of 2 A 
each. Each chunk produces an independent measure of the wavelength, PSF, and Doppler 
shift. The final measured velocity is the weighted mean of the velocities of the individual 
chunks. 
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Table [2] lists the complete set of 164 RVs for HD 1461, corrected for the solar system 
barycenter. The table lists the barycentric JD of observation center, RV, and internal 
uncertainty The internal uncertainty reflects only one term in the overall error budget, 
and results from a host of systematic errors from characterizing and determining the PSF, 
detector imperfections, optical aberrations, effects of under-sampling the iodine lines, 
etc. Two additional major sources of error are photon statistics and stellar jitter. The 
latter varies widely from star to star, and can be mitigated to some degree by selecting 
magnetically inactive older stars and by time-averaging over the star's unresolved low-degree 
surface p-modes. For most of the past 12.8 years, only single exposures at Keck were 
taken of HD 1461 at each epoch. Since these single exposures were much shorter than 
the characteristic time scale of low-degree surface p-modes on the star, they suffered from 
additional noise (stellar jitter). By 2008 July at Keck, we began p-mode averaging each 
observation, combining multiple shots of HD 1461 over a 5-10 minutes dwell at each epoch. 
All observations have been further binned on two-hour timescales. 

Although all of the observations used in this work were obtained with the HIRES at 
the Keck I telescope, 28 of the RVs in Table [2] were derived from publicly available spectra 
from the NASA Keck Observatory Archive. These velocities are marked with "Q01" in the 
observatory column in Table [2] The CCD format of the Q01 run was shifted by several 
angstroms relative to our long-term standard, introducing a zero-point velocity offset 
between Q01 and our standard set-up. This is an additional parameter to be determined in 
the fits discussed below. 



4. Photometry 



In addition to our RV observations, we acquired high-precision photometric 
observations of HD 1461 during 13 consecutive observing seasons from 1996 November to 
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2009 January with the T8 0.80 m automatic photometric telescope (APT) at the Fairborn 
Observatory. Our APTs can detect short-term, low-amplitude brightness variability in 
solar-type stars due to rotational modulation of the visibility of surface magnetic activity 
(spots and plages), as well as longer-term variations associated with the growth and decay 



of individual active regions and the occurrence of stellar magnetic cycles ((Henry. 



1999|). The 



photometric observations help us to establish w hether observed RV v ariations are caused 



stell ar activity or p l aneta ry reflex motion (e.g. 



and 



Paulson et al. 



Henry et al. 



2000ah . 



Queloz et al. 



(2001) 



(120041 ) have presented several examples of periodic RV variations in 



solar-type stars caused by photospheric spots and plages. The photo metric observations 



are also useful to search 



Charbonneau et al. 



2000). 



or transits of the planetary companions (e.g., 



Henry et al. 



2000b 



The T8 APT is equipped with a two-channel precision photometer that separates the 
Stromgren b and y passbands with a dichroic filter and takes simultaneous 30 s integrations 
with two electromagnetic interference (EMI) 9124QB bi-alkali photomultiplier tubes. The 
APT measures the difference in brightness between a program star and a nearby constant 
comparison star or stars. The typical precision of a single observation is approximately 
0.0015 mag, as measured for pairs of constant stars. The automatic tele scopes, photom eters, 



Henry 



(119991). 



observing procedures, and data reduction techniques are described in 
F urther details o n the d evelopment a n d ope ration of the automated telescopes can be found 
in 



Henry 


(1995a 


b) and 


Eaton et al. 


(2003) 



For HD 1461, we used HD 2361 (V = 7.89, B - V — 0.470, F2) as our primary 
comparison star. The individual Stromgren b and y differential magnitudes have been 
corrected for differential extinction with nightly extinction coefficients and transformed to 
the Stromgren system with yearly mean transformation coefficients. Since HD 1461 lies at 
a declination of —8°, the photometric observations from Fairborn were made at airmass 



- 7- 



1.3-1.8, which is somewhat higher than most observations. Therefore, to improve the 
precision of the measurements, we combined the Stromgren b and y differential magnitudes 
into a single (b + y)/2 passband. 

A total of 799 differential magnitudes from 13 observing seasons are plotted in the 
top panel of Figure [TJ The data scatter about their mean with a standard deviation of 
a = 0.0019 mag, which provides an upper limit to possible brightness variation in HD 1461. 
The data are plotted such that observations that are brighter/dimmer than the mean 
are indicated with positive/negative flux values. A Lomb-Scargle periodogram of the 
photometric measurements is shown in the bottom panel of Figure [Hand reveals a significant 
periodicity within the data at 444.5 days with an estim ated false alarm probability (FAP) 



of 3.3 x IP " 9 (cal culated using procedures described in 



Gilliland &: Baliunas 



Cumming 



(119871 ) and 



(J2004J)). This peak is near the tallest peak in the RV periodogram of the 
residuals of the circular two-planet fit discussed below. A second peak spans from 3000 
to 9000 days with an estimated FAP of 2.1 x 10 -6 . This broad peak is similar to the one 
present in the periodogram of the residuals of the (circular) one-planet fit also discussed 
below. The three horizontal lines in the bottom panel of Figure [1] and in all similar 
periodograms below represent FAPs of 0.1%, 1%, and 10% from top to bottom, respectively. 
We computed least-squares sine fits for the three RV planet candidates described below. 
Sinusoid fits with periods 5.77, 446.1, and 5017 days yield semi-amplitudes of only 0.00018, 
0.00070 and 0.00061 mag, respectively. We conclude that the lack of significant coherent 
photometric variability of HD 1461 supports planetary reflex motion as the cause of the 
5.77-day periodicity in the RV measurements. The weak periodicities observed at 445 
and 3000-9000 days, on the other hand, do warrant some caution in interpreting the RV 
periodicities at those periods as being due to planetary companions. 
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Fig. 1. — Top: differential photometry of HD 1461. Times are barycentric JD. Bottom: 
Lomb-Scargle periodogram of the photometry. 
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5. The Planetary System Orbiting HD 1461 

The RVs show a rms scatter of 3.80 ms 4 about the mean velocity. This significantly 
exceeds the combined scatter due to the underlying precision of both our measurement 
pipeline and the scatter expected in this star due to its predicted 1.59 ms" 1 level of stellar 
jitter. Figure [2] shows the RV data set from Table [2j An offset of -1.864 ms -1 (Q01 - Keck) 
has been applied between the two sets. Note that in general, observations prior to 2004 
August have larger internal uncertainties. The improvem ent in the inte r nal u ncertainties is 



Rivera et al. 



f|2005h . Also, the 



a result of the CCD upgrade discussed in some detail in 
observations prior to 2004 August appear to be offset from those taken after the CCD 
upgrade. This is an artifact of our sampling with low frequency a period that is not too 
distinct from one year during the first few years of our observations. As a result, during 
this time, we always sampled the period of HD 1461 c such that the star's RV is low. 

Figure [3] shows the periodogram of the RV data set (top panel). Power at each sampled 
period is proportional to the relative improvement (drop in xl) m the fit quality for a 
circular model versus a constant velocity model. The periodogram shows a number of 



significant signals, with the strongest peak occurring at a peri od of 5.77 d a ys 



Cumming 



he FAP 



(J2004J)) to be 



of this peak is estimated (adopting the procedure described in 
1.1 x 10~ 10 . Furthermore, this 5.77-day signal lies well away from the periods favored by 
the sampling window (Figure bottom panel), which produces spurious power at periods 
near 29.6, 186.3, 361.5, 147.8, and 82.4 days. Signals near any of these periods would be 
suspected of being artifacts of the observing scheduling. 

Based on the periodogram, with an assumed stellar mass of 1.022 M , we fit a planet 
of mass 7.4 M e and period 5.77 days on a circular orbit to the RV data. The presence 
of this planet (with a RV semi-amplitude of K = 2.60 ms -1 ) reduces the rms scatter of 
the velocity residuals to 3.43 ms -1 . Figure H] (top panel) shows the one-planet residuals 
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Fig. 2. — Relative RVs of HD 1461. Velocities based on spectra obtained by our group are 
in red, and those based on archived Keck spectra (Q01) are in blue. 
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Fig. 3. — Top: periodogram of the RV data set for HD 1461. The tallest peak is at 5.77 
days. Bottom: power spectral window of the RV data. 
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periodogram which has a strong peak at 4226 days. The 4226-day signal has a FAP of 
9.9 x 10 -16 and can be modeled with a companion with K — 2.77 ms" 1 and a mass of 0.22 
Mj up . The addition of this planet further reduces the rms scatter to 2.87 ms -1 . 

Figure H] (bottom panel) shows the two-planet residuals periodogram, assuming circular 
orbits, which has a peak at a period of 446.1 days. The FAP of the 446.1-day signal is 
4.5 x 10 -6 . If Keplerian, this periodicity can be ascribed to the presence of a 446.1-day, 
27.9 M e companion. This third component has K = 2.30 ms -1 . The three-planet model, 
assuming circular orbits, has an rms scatter of 2.41 ms 4 . 

Figure [5] shows the RV model for the three-planet fit with circular orbits with 
the observations overplotted. It is centered around the time of the high cadence Q01 
observations. 

Given the one-, two-, and three-planet models, we can look either for solutions in which 
the planetary orbits are circular or solutions where the eccentricities are allowed to float. 
Inclusion of eccentricities provides only a modest improvement to the orbital fits, leading 
us to conclude that a significant amount of additional Doppler velocity monitoring will be 
required to improve the eccentricity uncertainties. However, allowing the eccentricities to 
float for the two-planet fit dramatically reduces the significance of the 446-day periodicity 
observed in the two-planet residuals (see below). Additionally, if we fit for a second planet 
with period in the range ~390-450 days on an eccentric orbit, the significance of the 
long-period planet can also be reduced. In Tables [3] and HI we present our best-fit versions 
of the system under the assumption of circular orbits (Table [3]) and with the additional 
degrees of freedom provided by fully Keplerian trajectories (Table H]). For the orbital fits, 
we assume i = 90° and Q = 0°. The inclusion of planet-planet gravitational interactions 
in the fits were found to be unnecessary. Un certainties are base d on 1000 bootstrap trials 



following the procedure in Section 15.6 from 



Press et al. 



(Il992f ) . The standard deviations 
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Fig. 4. — Top: periodogram of the one-planet residuals, assuming a circular orbit, of the RV 
data set for HD 1461. Bottom: periodogram of the two-planet residuals, assuming circular 
orbits. 
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Fig. 5. — Zoomed-in view of the circular three-planet model for HD 1461 (black curve). 
Overplotted are the Keck observations in red and the Q01 observations in blue. 
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of the fitted parameters to the bootstrapped RV's were adopted as the uncertainties. 



The fitted mean anomalies are reported at epoch JD 2450366.936. T 
star is assumed to be 1.022 M^, the mean o f the isochrone masses of 



re mass of the host 



;i.018 ± 0.1 M Q ) and 



Sousa et al. 



Takeda et al. 



fl2008h 



(120071 ) ( 1.0261^ ). Our fitting was carried out 



with the publicly available Systemic Console (jMeschiari et al. 



2009] ) 



Figure [6] (top panel) shows the power spectrum of the velocity residuals for the 
two-planet fit with floating eccentricities. It is clear that the significance of the 446.1-day 
planet is dramatically reduced. 

Figure [6] (bottom panel) shows the power spectrum of the velocity residuals for the 
three-planet, circular fit. There is a peak near 18.3 days with a FAP of 4.7 x 10" 5 . Although 
using this period in a four-planet circular fit results in a significant improvement in xt> the 
rms decrease from 2.41 ms" 1 for the three-planet model to 2.28 ms -1 is not significant. 
Additionally, the fitted amplitude of 1.49 ms -1 is significantly smaller than the scatter 
around the model. 

In summary, the Keck HIRES RV data show strong evidence for at least two planets 
in orbit about the GOV star HD 1461. One of these has a short period of 5.77 days. The 
other has a large period of about 4000-6000 days, with a large uncertainty in both its 
period and eccentricity. There are also hints of a third planet in the system near 446 days, 
but the large uncertainties in the parameters of the second planet cast some doubt on the 
presence of a third planet in the system. Also, if we choose to fit for an eccentric planet 
with a period of ~390-450 days, we also find large uncertainties in its parameters. This 
would cast some doubt on the presence of the large period planet. The close match between 
the photometric and the Doppler periodicities near 445 days leads us to tread carefully in 
ascribing a planetary origin to our 446-day RV signal. We note, however, that a 446-day 
period is too long to be associated with the stellar rotation, given the star's measurable 
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Fig. 6. — Top: periodogram of the two-planet residuals, with floating eccentricities, of the 
RV data set for HD 1461. Bottom: periodogram of the three-planet residuals, assuming 
circular orbits. 



-17- 



rotational velocity component, Vsini = 1.6 km s -1 . If the 446-day Doppler signal does arise 
from stellar activity, then the periodicity would need to be primarily associated with the 
lifetimes of the active regions rather than with modulation induced by stellar rotation. 
Further observations will be needed to refine the parameters of the second planet before a 
definitive detection claim can be made for a third planet. That said, in Figure [7| we show 
the barycentric reflex velocity of the host star due to each of three individual companions 
on circular orbits in the system. In each panel, the velocities are folded at the period of 
each corresponding planet. 



6. Dynamical Analysis 



It is useful to verify, via numerical integration, whether the planetary configurations 
listed in Tables |3] and H] are dynamically stable. Such an analysis is particularly useful 
in giving rough bounds on the allowed coplanar inclinations relative to the line of sight 
to Earth. For simulations that do not include the effects of tidal dissipation, Newtonian 
parameters were use d for the initial s tates for long-term integrations. The MERCURY 



19991 ) was used for the simulations with a time step of 0.1 



integration package (iChambers 
day. The first orde r post -Newtonian term in the star's potential was also included, as in 



Lissauer fc Rivera 



(j200lh . 



If the three orbits are assumed to be initially circular, with periods, masses and 
mean anomalies given in Table [3j the system is stable for at least 20 Myr. Additionally, 
assuming the system to be coplanar, if we set the inclination to the sky plane to various 
values from i = 90° all the way down to i = 1° and perform a Newtonian fit for the other 
nine parameters (three parameters per planet plus the two offsets), xl does not change 
significantly from the nominal i = 90° fit. We also find the i — 1° fit to result in a system 
that is stable for at least 20 Myr. For this inclination, the fitted masses exceed 1.4, 5.1, 
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Fig. 7. — Top: radial velocity of HD 1461 due to planet b folded at 5.77 days. Center: radial 
velocity of HD 1461 due to planet c folded at 446.1 days. Bottom: radial velocity of HD 1461 
due to planet d folded at 5017 days. In each panel, the effect of the other two planets has 
been subtracted out. The curves represent the model velocities due to each respective planet. 
The Keck observations are shown in red, and the Q01 velocities are shown in blue. 
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and 16 Mj up . This system is stable because of the small eccentricities. Thus, under the 
assumption that the system is coplanar and the orbits are (nearly) circular, we cannot place 
a lower bound on the inclination of the system. 

The parameters of the floating-eccentricity version of the HD 1461 system given 
in Table H] were also used as the initial input conditions for a 10 7 yr simulation. The 
three-planet configuration is disrupted in less than 1.4 Myr. 

If additional RV measurements point to secure non-zero eccentricities for the HD 1461 
planets, then studies of the long-term dynamical evolution of the system should take the 
possibility of tidal dissipation in planet b into account. 



7. Check for Transits by Companion b 

We performed a simple search for transits of planet b in the HD 1461 photometry. The 
RV's determine the period of the inner planet to be 5.77 day. The corresponding mean 
anomaly of the planet, likewise, is M = 79.045°. The stellar radius is taken to be 1.1 R Q , 
the mean of the values in Table [TJ 

If HD 1461 b's orbital plane allows for transits, then sin % ~ 1, and hence Mf, = 7.4 M®. 
If we assume that the planet migrated inw ard from beyond the ice line, its composition is 



likely dominated by water. The models of iFortney et al. 



(120071 ) suggest an R = 2.7 R Q 



radius for such a planet, leading to a central transit depth of d ~ 0.05%, or d ~ 0.5 mmag. 
If the planet has a massive atmosphere, its transit depth will be larger still. 

At the time of the first photometric data point, T pl = JD2450393.7339, our orbital 
model indicates a mean anomaly M = 310.255° for planet b. For a circular orbit, transits 
occur at M — 90°. This means that relative to phase <p — 0.0 referenced to T p \, transits 
are centered at <ft = 0.388. The planet's orbital velocity is i>& ~ 120 kms -1 which, for our 
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R = 1.1 Rq star, implies a transit duration of r = 12000 sec = 3.4 hours. Assuming central 
transits, planet b would be observable in transit 2.4% of the time. 

We assumed a simple "top-hat" step- function model for the transit, with the depth, 
period, and phase given above. The photometric data set contains 799 individual 
measurements, with a = 0.00188 mag. Phase coverage at the 5.77-day period of planet b 
is excellent. One thus expects N = 19 ± 4 points to lie in the transit window, indicating 
that, given the data set and the presence of a transiting planet b, one can expect to 
detect a transit with signal-to-noise ratio, S/N~l. We therefore conclude that the present 
photometric data set is insufficient to make a definitive call as to whether transits are 
occurring. Follow-up with high-precision, high-cadence photometry from either ground or 
space is therefore warranted. 



8. Discussion 

We present evidence for at least two (and possibly three) low-mass planets orbiting 
the nearby star HD 1461. With a RV semi-amplitude K = 2.8 ± 0.3 ms _1 , the 
M smi = 7.8 ± 0.8 M® inner planet HD 1461b is among the very lowest-amplitude 
companions yet detected using the Doppler velocity technique. 

The HD 1461 syste m is thus another ne arby case that joins the emerging population 



of planets postulated by 



Mayor et al. 



(120091 ). who inferred that about 30% of solar-type 



stars in the immediate galactic neighborhood are accompanied by Neptune (or lower) mass 
planets having orbital periods of 50 days or less. 

With a period of only 5.77 days, HD 1461b has a non-negligible P ~ 8% probability 
of transiting its parent star. While our phase-folded ground-based photometry does not 
have the requisite cadence and S/N to detect such transits, it would be readily possible to 
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determine whether transits occur by making high-cadence, high-precision observations that 
span the transit window (see, e.g. Johnson et al. 2009). HD 1461 b's mass is very similar 
to that of CoRoT-7b, which has recently been determined to have a density similar to that 
expected for a rocky planet (Queloz et al. 2009). It would be of great interest to learn 
whether HD 1461 b is similarly dense, or whether its composition is more reminiscent of 
ice-giant planets such as Neptune, Gliese 436b, and HAT-P-llb. Such a transit could be 
detected from space using, for example, the Warm Spitzer platform. 

As RV data bases grow in the monitoring of chromospherically quiet, nearby stars, 
systems like HD 1461, are becoming increasingly common. With each new system, the 
evidence is growing stronger that Super-Earths and other low-mass planets are common 
around nearby Sun-like and cooler stars. It is thus only a matter of time and adequate 
cadence before Super-Earth planets are found in the habitable zones of nearby stars. This 
is only a first reconnaissance of this fascinating and quite nearby system. As RV data bases 
grow for this star, the orbital ephemerides of these planets will become better determined, 
and more planets will probably be revealed. 
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Table 1. Stellar Parameters for HD 1461 



Parameter 



Value 



Reference 



Spec. Type 
Mass (Mo) 



Radius (Rq) 

Luminosity (1/q) 

Distance (pc) 
V sini (kms -1 ) 
logi^ K 



P ro t (days) 
age (Gyr) 



[Fe/H] 

T ctt (K) 
logs 



GOV 
1.08 ± 0.04 
1 026+ 040 

l.UZD_g 03 q 

1.018 ± 0.1 
1.095 ± 0.026 
1.11 ± 0.04 
1.197±0.113 
1.188 ±0.017 
23.4 ±0.5 
1.6 
-5.03 
-5.04 
-5.00 
-5.00 
29 
6.3 

4.2 H 



>+l.7 
-2.1 



7 i o +1.40 
'■ iz -1.56 

0.18 
0.20 + 0.01 

0.19 + 0.01 
5765 ± 18 
5808 

4.37 + 0.03 

4.39 

4.38 + 0.03 



Cenarro et al. (2007) 

Valenti fc Fischer (2005) 

Takeda et al. (2007) 

Sousa et al. (2008) 

Valenti fc Fischer (2005) 

Takeda et al. (2007) 

Valenti fc Fischer (2005) 

Sousa et al. (2008) 

F^CTr^rnan_et_aL_ (1997) 

Valenti fc Fischer (2005) 

Wright et al. (2004) 
Hall et al. 



(2007) 
(2009) 



Hall etaj. 

This work 
Wright et al. (2004) 
Wright et al. 



Valenti . 



(2004) 

Valenti fc Fischer (2005) 

Takeda et al. (2007) 

Valenti fc Fischer (2005) 

Cenarro et al. (2007) 

Sousa et al. (2008) 

Fischer (2005) : Sousa et al. (2008) 

Cenarro et al. (2007) 



Valenti fc Fischer (2005) ; Takeda et al. 

Cenarro et al. 
Sousa et al. 



(2007) 
(2007) 
(2008) 



Table 2. Radial Velocities for HD 1461 



Barycentric JD RV Error Observatory 

(-2450000) (ms- 1 ) (ms- 1 ) 



366.93641 


-0.21 


1.36 


K 


715.00792 


-9.17 


1.11 


K 


785.76804 


-14.14 


1.22 


K 


1010.09403 


-5.98 


1.17 


K 


1068.92515 


-5.81 


1.14 


K 


1173.72568 


-0.48 


1.10 


K 


1343.08899 


-2.04 


1.12 


K 


1374.10891 


-4.86 


1.27 


K 


1438.88093 


-3.75 


1.55 


K 


1550.73862 


-4.91 


1.25 


K 


1883.80712 


-0.32 


1.34 


K 


2098.12817 


-7.96 


1.51 


K 


2236.72882 


-0.08 


1.30 


K 


2489.05617 


-2.67 


1.51 


K 


2537.92834 


-3.48 


1.57 


K 


2575.84461 


-1.98 


1.54 


K 


2829.05914 


-7.47 


1.43 


K 


2899.03253 


-2.51 


1.37 


K 


3196.06321 


2.86 


1.28 


K 


3238.01691 


0.50 


0.21 


K 


3238.90057 


2.88 


0.42 


K 


3239.94328 


6.80 


0.44 


K 


3240.99951 


6.49 


0.49 


K 


3241.11480 


7.06 


0.65 


K 


3301.85378 


-4.94 


0.65 


K 


3302.81034 


-0.65 


0.61 


K 


3338.75610 


1.71 


1.08 


K 


3339.75947 


-1.18 


0.66 


K 


3369.75405 


-6.98 


0.64 


K 


3397.70379 


3.24 


0.58 


K 


3398.73620 


-0.30 


0.61 


K 
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Table 2 — Continued 



Barycentric JD RV Error Observatory 

(-2450000) (ms- 1 ) (ms^ 1 ) 



3399.69869 


-2.37 


0.64 


K 


3400.70946 


4.84 


0.64 


K 


3547.12033 


3.03 


0.63 


K 


3548.10688 


3.59 


0.96 


K 


3549.13414 


0.06 


1.14 


K 


3550.11494 


3.57 


0.61 


K 


3551.12904 


4.00 


0.62 


K 


3552.12582 


4.68 


0.58 


K 


3603.04546 


7.42 


0.63 


K 


3604.06836 


10.32 


0.75 


K 


3605.04697 


6.44 


0.64 


K 


3692.92330 


-2.41 


0.61 


K 


3693.84499 


1.79 


0.54 


K 


3694.78313 


1.37 


0.58 


K 


3695.70050 


4.03 


0.53 


K 


3695.82845 


4.85 


0.59 


K 


3696.76371 


5.63 


0.57 


K 


3723.82034 


-1.81 


0.60 


K 


3724.74409 


1.77 


0.63 


K 


3746.73490 


-1.19 


0.69 


K 


3747.79465 


-0.46 


0.97 


K 


3748.72482 


1.91 


0.68 


K 


3749.72744 


-0.78 


0.62 


K 


3775.71584 


-2.90 


0.67 


K 


3776.70867 


2.10 


0.65 


K 


3777.71251 


1.53 


0.80 


K 


3778.71106 


-2.19 


0.86 


K 


3779.73772 


-2.38 


0.81 


K 


3927.10455 


0.42 


0.62 


K 


3928.05300 


1.31 


0.85 


K 


3959.10907 


-1.33 


0.62 


K 
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Table 2 — Continued 



Barycentric JD RV Error Observatory 

(-2450000) (ms- 1 ) (ms^ 1 ) 



3961.04337 


0.02 


0.64 


K 


3962.01510 


2.99 


0.61 


K 


3981.90046 


-3.83 


0.65 


K 


3982.98884 


2.04 


0.66 


K 


3983.87130 


2.26 


0.58 


K 


3984.97007 


6.59 


0.65 


K 


4083.81187 


2.61 


0.73 


K 


4084.77778 


3.39 


0.73 


K 


4129.74967 


8.16 


0.66 


K 


4131.71246 


4.76 


0.75 


K 


4279.11113 


0.46 


1.19 


K 


4280.10709 


-1.69 


1.26 


K 


4286.10589 


-5.34 


1.22 


K 


4295.08454 


-2.01 


1.15 


K 


4305.06828 


-5.52 


0.61 


K 


4306.02305 


-1.97 


0.83 


K 


4306.99662 


-3.53 


0.94 


K 


4307.10265 


-0.77 


0.83 


K 


4308.07916 


-0.27 


0.71 


K 


4309.07893 


0.58 


0.85 


K 


4310.02797 


-4.99 


0.80 


K 


4310.12123 


-2.37 


1.19 


K 


4311.00946 


-6.55 


0.84 


K 


4311.10912 


-2.30 


1.23 


K 


4312.00404 


-4.29 


0.82 


K 


4312.11147 


-2.84 


0.82 


K 


4313.00139 


-3.43 


0.87 


K 


4313.10575 


-2.31 


0.86 


K 


4313.99949 


-0.07 


0.85 


K 


4314.10977 


0.60 


0.90 


K 


4315.12175 


-1.56 


0.86 


K 



-30- 



Table 2 — Continued 



Barycentric JD RV Error Observatory 

(-2450000) (ms- 1 ) (ms^ 1 ) 



4319.01064 


3.70 


0.70 


K 


4336.05526 


1.63 


0.75 


K 


4337.10146 


4.33 


0.81 


K 


4343.93370 


5.91 


1.21 


K 


4396.76696 


-2.48 


1.18 


K 


4397.82001 


-4.10 


1.23 


K 


4398.84808 


-2.28 


1.12 


K 


4427.86323 


0.36 


0.67 


K 


4429.77428 


2.62 


1.30 


K 


4430.77076 


2.35 


1.24 


K 


4454.82508 


-2.21 


1.20 


K 


4456.80064 


-2.24 


1.25 


K 


4460.78595 


-4.89 


1.35 


K 


4461.79834 


-6.06 


1.13 


K 


4464.80653 


4.89 


1.37 


K 


4492.70875 


5.97 


1.29 


K 


4634.11567 


-6.90 


0.92 


K 


4635.07070 


-5.12 


0.88 


K 


4636.09364 


-3.98 


0.87 


K 


4637.10870 


-6.08 


0.88 


K 


4638.11047 


-2.74 


0.91 


K 


4639.09420 


-4.49 


0.88 


K 


4640.10790 


-4.80 


0.84 


K 


4641.10062 


-7.37 


0.86 


K 


4642.07466 


1.06 


0.94 


K 


4644.12501 


1.89 


0.89 


K 


4717.94972 


-2.96 


0.83 


K 


4719.00539 


0.38 


0.73 


K 


4720.00787 


-4.48 


0.66 


K 


4720.96693 


-4.91 


0.77 


K 


4721.98581 


-5.73 


0.79 


K 



Table 2 — Continued 



Barycentric JD RV Error Observatory 

(-2450000) (ms- 1 ) (ms^ 1 ) 



4722.89844 


-0.54 


0.71 


K 


4723.96804 


0.62 


0.73 


K 


4724.96948 


-0.12 


0.82 


K 


5022.12042 


-5.09 


0.52 


K 


5023.08514 


-3.59 


0.85 


K 


5024.11008 


-1.22 


0.49 


K 


5025.12159 


-4.13 


0.58 


K 


5050.06798 


-2.46 


0.49 


K 


5051.09955 


-2.43 


0.46 


K 


5052.04858 


1.40 


0.55 


K 


5053.08104 


-0.45 


0.56 


K 


5054.04263 


4.05 


0.59 


K 


3713.70121 


0.45 


0.70 


Q01 


3713.76642 


1.42 


0.59 


Q01 


3713.82984 


3.44 


0.60 


Q01 


3714.69699 


0.24 


0.61 


Q01 


3714.75572 


0.56 


0.58 


Q01 


3714.82460 


0.22 


0.56 


Q01 


3715.68459 


-1.65 


0.66 


Q01 


3715.72791 


0.10 


0.70 


Q01 


3715.77573 


-2.24 


0.64 


Q01 


3715.82808 


-0.68 


0.63 


Q01 


3716.72884 


-3.10 


0.59 


Q01 


3716.77010 


-1.49 


0.63 


Q01 


3716.82182 


-1.79 


0.72 


Q01 


3717.79314 


-4.50 


0.56 


Q01 


3726.72138 


-1.00 


0.61 


Q01 


3726.77494 


2.24 


0.59 


Q01 


3727.71647 


-1.39 


0.54 


Q01 


3727.77202 


-4.18 


0.56 


Q01 


3727.81049 


-2.38 


0.58 


Q01 
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Table 2 — Continued 



Barycentric JD 


RV 


Error 


Observatory 


(-2450000) 


(ms- 1 ) 


(ms- 1 ) 




3728.73194 


-5.62 


0.54 


Q01 


3728.79946 


-2.87 


0.59 


Q01 


3728.84615 


-2.39 


0.58 


Q01 


3728.84971 


0.18 


0.60 


Q01 


3729.70727 


1.63 


0.85 


Q01 


3729.81947 


0.00 


0.60 


Q01 


3730.71913 


2.23 


0.60 


Q01 


3730.77245 


2.82 


0.56 


Q01 


3730.80874 


2.49 


0.54 


Q01 



Table 3. Circular Solutions (Epoch JD 2450366.936) 



Planet 


Period 


K e 




M 


Msini 


a 




(d) 


(m s -1 ) 


(deg) 


(deg) 




(AU) 








1-planct; \i = 


19.249; 


rms=3.43ms" 


-l 




HD1461 b 


5.7726±0.0026 


2.6±0.4 0.0 


n/a 


79±60 


7.4±1.2 


0.063437±0.000019 






2-plancts; \t = 


13.151; 


rms=2.87ms 


-l 




HD1461 b 


5.7718±0.0010 


2.8±0.3 0.0 


n/a 


44±36 


7.8±0.8 


0.063431±0. 000008 


HD1461 d 


4019±1433 


2.8±1.9 0.0 


n/a 


285±47 


69.9±78.9 


4.98±0.99 






3-plancts; \i = 


10.093; 


rms=2.41 m s 


-l 




HD1461 b 


5.7718±0.0010 


2.7±0.2 0.0 


n/a 


45±37 


7.6±0.7 


0.063431±0.000008 


HD1461 c 


446±9 


2.3±0.4 0.0 


n/a 


82±52 


27.9±4.9 


1.151±0.016 


HD1461 d 


5017±1171 


3.2±1.4 0.0 


n/a 


319±36 


87.1±51.8 


5.78±0.82 
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Table 4. Eccentric Solutions (Epoch JD 2450366.936) 



Planet 


Period 
(d) 


K e u> 
(ms -1 ) (deg) 


M 
(deg) 


M sin i 


a 

(AU) 








1-planet; \l = 19.349; rms= 


3.43 ms- 1 






HD1461 b 


5.7727±0.0025 


2.7±0.6 0.14±0.19 58±69 


26±71 


7.6±1.2 


0.063438±0.000018 






2-planets; xt = 12.768; rms= 


:2.64 ms- 1 






HD1461 b 
HD1461 d 


5.7720±0.0011 
7000±200000 


2.8±0.4 0.06±0.14 28±70 
4±24 0.51±0.24 319±32 


24±67 
210±83 


8.1±0.9 
101±357 


0.063432±0.000008 
7±42 






3-planets; \1 = 8.983; rms= 


2.24ms- 1 







HD1461 b 
HD1461 c 
HD1461 d 



5.7722±0.0011 
454±4 
5000±90000 



2.8±0.3 
2.8±2.7 
4±70 



0.04±0.01 186±63 234±77 8.1±0.7 0.063434±0.000008 
0.74±0.13 87±28 237±28 22.9±9.8 1.165±0.008 
0.16±0.29 326±79 136±85 97±1161 5±18 



